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^ ! ABSTRACT 

' Aims. To investigate the first high resolution optical spectrum of the B-type star, LS III -h52°24, 

O ' 

identified as the optical counterpart of the hot post-AGB candidate IRAS 22023+5249 (122023). 
■ Methods. We carried out detailed identifications of the observed absorption and emission fea- 

. tures in the high resolution spectrum (4290 - 9015 A ) of 122023 obtained with the Utrecht 



Echelle Spectrograph on the 4.2m William Herschel Telescope. Using Kurucz's WIDTH9 pro- 
gram and the spectrum synthesis code, SYNSPEC, we determined the atmospheric parameters 
and abundances. The photospheric abundances were derived under the LTE approximation. The 



' NEBULAR package under IRAF was used to estimate the electron temperature (Te) and the elec- 

tron density (N^) from the [N II] and [S II] lines. 

Results. We estimated Teff =24000 K, log g=3.0, f,=7 kms"'. The derived CNO abundances 
suggest an evolved star with C/O < 1. P-Cygni profiles of hydrogen and helium indicate on- 
going post-AGB mass loss. The presence of [N II] and [S II] lines and the non-detection of 
[O III] indicate that photoionisation has just started. The derived nebular parameters (Tc=7000 K, 
Nc= 1 .2x lO'' cm"') are comparable to those measured in young compact planetary nebulae (PNe). 
The nebular expansion velocity was estimated to be 17.5 kms"'. 

Conclusions. The observed spectral features, large heliocentric radial velocity (-148.31 ± 0.60 
kms"'), atmospheric parameters and chemical composition indicate that 122023, at a distance 
of 1.95 kpc, is an evolved post-AGB star belonging to the old disk population. The nebular 
parameters suggest that the central star may be evolving into a compact, young PN, similar to 
Hen3-I357. 



Key words. Stars: AGB and post-AGB - Stars: early-type - Stars: abundances - Stars: evolu- 
tion 
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1. Introduction 



The discovery of high galactic latitude cool (F, G, K) and hot (OB) post-AGB supergiants, e.g. 
HD 161796 (Parthasarathy & Pottasch 1986) and LSI! +34°26 (Parthasarathy 1993) indicated that 
K, G, F, A, OB post-AGB supergiants form an evolutionary sequence in the transition region 
from the tip of the AGB into early stages of planetary nebulae (PNe). Since then, several cool 
and hot post-AGB candidates have been identified (Pottasch & Parthasarathy 1988, Parthasarathy 
& Pottasch 1989, Garci'a-Lario et al. 1997a; Parthasarathy et al. 2000a). Gauba & Parthasarathy 
(2003, 2004) analysed the UV spectra and circumstellar dust envelopes of several hot post-AGB 
stars from the above lists, including LS 111 +52°24. High resolution optical spectra of only a few 
hot post-AGB stars have been studied (Sarkar et al. 2005 and references therein). To ascertain their 
evolutionary nature, atmospheric parameters and chemical compositions, there is a clear need to 
study a bigger sample of these objects. 

The optically bright B-type star, LS HI +52°24 identified with the IR source 
IRAS 22023+5249 (hereafter, 122023) has far-IR colors similar to PNe (Table 1). It is listed in 
WackerUng's (1970) Catalog of early-type emission-line stars. In this paper, we present the first 
high resolution optical spectrum of 122023. Recent, ground-based, high resolution, near-IR im- 
ages, have shown H2 emission arising close to the central star, possibly in a bipolar morphology 
(Volk et al. 2004). Kelly & Hrivnak (2005) also detected Hydrogen Br^ 2.166 fxm line in emission. 
Gauba & Parthasarathy (2004) indicated the presence of weak amorphous (10.8/im) and crystalline 
(33.6/zm) silicate features in the ISO spectrum of the object. However, another look at the ISO 
spectrum suggests that this identification was possibly in error. The ISO spectrum of 122023 is 
very noisy and displays only a featureless continuum with no 10.8 jum and 33.6 jum features. 

Table 1. Details of the star 



IRAS 


Name 


RA 


DEC 


1 


b 


Sp. Type 


V 


B-V 




IRAS Fluxes (Jy) 








2()()() 


2000 






Optical 






12 /.mi 


25 //m 


60 /mi 


100 jum 


22023+5249 


LSm +52 24 


22:04:12.30 


+53:04:01.4 


99.30 


-1.96 


B" 


12.52" 


0.69" 


1.02 


24.69 


14.52 


3.93L 



"Spectral type is from the SIMBAD database. 'Photometry is from Hog et al. (2000). L flag indicates that the quoted IRAS flux 



density is an upper limit. 



2. Observations 

122023 was observed on 14 July 2001 using the Utrecht Echelle Spectrograph (UES) on the 4.2m 
Wilham Herschel Telescope (WHT) at the Roque de los Muchachos Observatory on La Palma. 
Observations were made with the 31.6 lines/nrai echelle grating (E31), SITel CCD (2048 x 2048 
pixels of 24 //m), a sht width of l"on the sky and central wavelength of 5 5 00 A resulting in a 
resolving power of R ~ 52000. The wavelength coverage was 4290-4735A , 4760-5553A and 
5607-9015A . A Th-Ar comparison lamp was used for wavelength cahbration. 

Send offprint requests to: G. Sarkar 
* Based on observations obtained with the 4.2m William Herschel Telescope and the Utrecht Echelle 

Spectrograph of the Roque de Los Muchachos Observatory, La Palma, Spain 
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3. Analysis 

The spectrum was extracted using standard reduction procedures in the IRAF package. The reduc- 
tion steps included overscan correction, trimming, bias subtraction, flat field correction, coiTection 
for scattered light and wavelength calibration. The reduced spectrum was continuum-normalised. 
The final S/N varied from 30 in the blue to 60 towards the red end of the spectrum. 

Equivalent widths (W^) of the absorption and emission lines were measured. Deblending 
was done whenever required to obtain Gaussian fits to the blended line profiles. The complete 
continuum-normalised spectrum of 122023 (Fig. 1) is presented in the appendix. This spectrum 
would be useful for future observers since post-AGB stars show both short and long-term vari- 
ability in the absorption and emission line strengths and profiles. P-Cygni profiles detected in these 
stars are also expected to vary as the stellar wind and post-AGB mass loss rates may show vari- 
ation as the star evolves. The line identifications (Tables 2-5) are based on the Moore multiplet 
table (1945) and flie lineUsts of Parfliasarathy et al. (2000b), Klochkova et al. (2002) and Sarkar et 
al. (2005). Unidentified fines are denoted by "UN". Night sky emission lines denoted by "atmos." 
were identified from Osterbrock et al. (1996). The laboratory wavelengths, log (gf) values and ex- 
citation potentials (x) have been extracted from the Kurucz (CD-ROM 23) linelist GFALL. Ivan 
Hubeny and Thierry Lanz have compiled the Kurucz linelists with improved oscillator strengths 
from the NIST Atomic Spectra Database. For wavelengths below 7500 A we have used their data 
which may be retrieved from http://nova.astro.umd.edu/Synspec43/synspec-frames-data.html, 

3. 1 . Description of tiie spectrum 

The absorption and emission lines in the spectrum of 122023 are similar to those detected in other 
hot post-AGB stars (Parthasarathy et al. 2000b, Klochkova et al. 2002 and Sarkar et al. 2005). 
In addition to the O I triplet, absorption lines of He I, C II, C III, N II, O II, Ne I, Mg II, Al III, 
Si III were identified. Emission lines of C II, N II, O I, [O I], Al III, Si II, Fe I, Fe II, [Fe II] and 
Fe III were also identified in the spectrum of the star. The Nal D 1 and D2 lines show a complex 
structure. The presence of low excitation nebular lines of [N II] and [S II] and the absence of [O III] 
5007A suggest that photoionization has just started, shock excitation may also be playing a role. 
Balmer lines of Ha, and Hy show P-Cygni profiles indicating ongoing post-AGB mass loss. 
Some He I, C II and Fe III lines were also found to have P-Cygni profiles. 

3.2. Reddening 

For the derived effective temperature (Tgff) and gravity of the star (Sec. 3.8), we may assume a 
BlIII spectral type with intrinsic B-V ((B-V)o) = -0.26 (Schmidt-Kaler 1982). The observed 
B-V for 122023 is 0.69 (Hog et al. 2000). The total E(B-V) towards the star is therefore estimated 
to be 0.95. At the low galactic latitude (b - -1.96) of the star, the Diffuse Infrared Background 
Experiment (DIRBE)/IRAS dust maps (Schlegel et al. 1998) do not give reliable estimates of the 
interstellar extinction. The interstellar E(B-V) was estimated to be 0.63 based on the strength of 
the DIB at 5780.410A (Sec. 3.6). This suggests that 122023 suffers from significant circumstellar 
reddening. The UV spectrum of the star obtained with the International Ultraviolet Explorer (lUE) 
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Table 2. Absorption Unes in IRAS 22023+5249 



'lobs. 


'llab. 


Ident. 




log (gf) 


X 






(A) 


(A) 




(A) 




(eV) 


(A) 


km s"' 


4314.874 


4317.139 


on(2) 


0.1609 


-0.386 


22.95-25.82 


2.265 


-142.53 


4317.277 


4319.630 


on(2) 


0.2197 


-0.380 


22.96-25.83 


2.353 


-148.55 


4343.245 


4345.560 


011(2) 


0.1448 


-0.346 


22.96-25.81 


2.315 


-144.95 


4345.141 


4347.413 


011(16) 


0.0480 


0.024 


25.64-28.49 


2.272 


-141.91 


4347.137 


4349.426 


011(2) 


0.2489 


0.060 


22.98-25.83 


2.289 


-143.01 


4348.903 


4351.260 


OII(16) 


0.0483 


0.227 


25.64-28.49 


2.357 


-147.64 


4364.653 


4366.895 


on(2) 


0.1130 


-0.348 


22.98-25.82 


2.242 


-139.15 


4385.301 


4387.929 


Hel(51) 


0.3722 


-0.883 


21.20-24.03 


2.628 


-164.80 


4412.525 


4414.899" 


011(5) 




0.172 


23.42-26.23 




blend 


4414.702 


4416.975* 


011(5) 


0.1319 


-0.077 


23.40-26.21 


2.384 


blend 


4435.158 


4437.551 


Hel(50) 


0.0786 


-2.034 


21.20-23.99 


2.393 


-146.91 


4478.638 


4479.885 


A1II1(8) 


0.0524 


0.900 


20.77-23.53 








+ 4479.971 


Aini(8) 




1.020 


20.77-23.53 








+ 4481.126 


Mgn(4) 




0.740 


8.86-11.62 


2.488 


-151.70 


4550.109 


4552.622 


Silll(2) 


0.3600 


0.181 


19.00-21.72 


2.513 


-150.73 


4565.371 


4567.840 


Silll(2) 


0.3225 


-0.039 


19.00-21.72 


2.469 


-147.28 


4572.278 


4574.757 


Silll(2) 


0.1638 


-0.509 


19.00-21.71 


2.479 


-147.70 


4588.513 


4590.974 


011(15) 


0.1206 


0.350 


25.64-28.34 


2.461 


-145.95 


4593.736 


4596.177 


011(15) 


0.1507 


0.200 


25.64-28.34 


2.441 


-144.46 


4627.829 


4630.539 


Nn(5) 


0.0983 


0.094 


18.47-21.14 


2.710 


-160.70 


4636.365 


4638.856 


011(1) 


0.1682 


-0.332 


22.95-25.62 


2.491 


-146.23 


4639.279 


4641.810 


011(1) 


0.3179 


0.055 


22.96-25.63 


2.531 


-148.71 


4644.997 


4647.418 


Clll(l) 


0.1181 


0.070 


29.51-32.18 


2.421 


-141.41 


4646.569 


4649.135 


on(i) 


0.3387 


0.308 


22.98-25.65 


2.566 


-150.71 


4648.297 


4650.16 


cni(i) 


0.2709 








blend 




+ 4650.841 


on(i) 














+ 4651.35 


Clll(l) 












4659.126 


4661.632 


011(1) 


0.2431 


-0.278 


22.96-25.62 


2.506 


-146.40 


4671.265 


4673.733 


011(1) 


0.0644 


-1.090 


22.96-25.61 


2.468 


-143.55 


4673.706 


4676.235 


on(i) 


0.1598 


-0.394 


22.98-25.63 


2.529 


-147.38 


4693.804 


4696.353 


on(i) 


0.0346 


-1.380 


22.98-25.62 


2.558 


-148.53 


4696.824 


4699.218 


011(25) 


0.0944 


0.270 


26.21-28.84 


2.394 


-137.96 


4702.873 


4705.346 


OIl(25) 


0.0621 


0.477 


26.23-28.86 


2.473 


-142.80 


4817.027 


4819.712 


Silll(9) 


0.0556 


0.750 


25.96-28.53 


2.685 


-152.26 


4826.346 


4828.951 


Silll(9) 


0.0344 


1.090 


25.97-28.53 


2.605 


-146.97 


4904.197 


4906.830"^ 


011(28) 




-0.161 


26.29-28.81 


2.633 


-146.11 


4918.616 


4920.35 


Hel(49) 


0.4414 




21.13-23.64 




blend 




+ 4921.931 


Hel(48) 




-0.435 


21.20-23.72 






4921.878 


4924.529 


OIl(28) 


0.0851 


0.074 


26.29-28.80 


2.651 


-146.63 



"011(5) 4414.899A is blended with FeII(32) A emission. *0n(5) 4416.975A is blended with the absorption component of 
Fein(114) 449.596 P-Cygni profile. ^On(28) 4906.88 A is blended with [FeII](20F) 4905.35 A emission. 
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Table 2. Absorption Unes in IRAS 22023+5249 



'lobs. 


'llab. 


Ident. 




log (gf) 


X 






(A) 


(A) 




(A) 




(eV) 


(A) 


km s"' 


4940.268 


4943.005 


011(33) 


0.0825 


0.239 


26.54-29.05 


2.737 


-151.24 


5000.141 


5002.703 


Nn(4) 


0.0328 


-1.022 


18.45-20.92 


2.562 


-138.77 


5002.414 


5005.150 


NII(19) 


0.0336 


0.594 


20.65-23.13 


2.736 


-149.12 


5007.895 


5010.621 


NII(4) 


0.0851 


-0.607 


18.45-20.92 


2.726 


-148.34 


5042.253 


5044.8 


CII(35) 


0.1148 








blend 




+ 5045.098 


NII(4) 












5044.65 


5047.2 


CII(35) 


0.0825 








blend 




+ 5047.736 


Hel(47) 












5130.277 


5132.947 


CII(16) 


0.0376 


-0.211 


20.69-23.10 




blend 




+ 5133.282 


CII(16) 




-0.178 


20.69-23.10 






5140.722 


5143.495 


CII(16) 


0.0640 


-0.212 


20.69-23.10 




blend 


5142.407 


5145.165 


CII(16) 


0.0523 


0.189 


20.70-23.10 




blend 


5148.32 


5151.085 


CII(16) 


0.0446 


-0.179 


20.70-23.10 


2.765 


-146.16 


5216.486 




UN 


0.0471 










5486.97 


5487.43 


DIB 


0.0343 






0.460 


-10.28 


5663.356 


5666.630 


NII(3) 


0.1400 


-0.045 


18.45-20.64 


3.274 


-158.46 


5672.88 


5676.02 


NII(3) 


0.1433 


-0.367 


18.45-20.63 


3.14 


-151.09 


5676.262 


5679.56 


Nn(3) 


0.2223 


0.250 


18.47-20.65 


3.298 


-159.33 


5682.922 


5686.21 


Nn(3) 


0.0382 


-0.549 


18.45-20.63 


3.288 


-158.60 


5693.475 


5695.920 


cni(2) 


0.1418 


0.017 


32.08-34.26 








+ 5696.604 


AIIII(2) 




0.230 


15.63-17.81 


3.129 


-149.91 


5707.504 


5710.770 


NII(3) 


0.0552 


-0.518 


18.47-20.64 


3.266 


-156.70 


5719.637 


5722.730 


A1III(2) 


0.0786 


-0.070 


15.63-17.80 


3.093 


-147.27 


5736.688 


5739.734 


Sini(4) 


0.2575 


-0.160 


19.71-21.87 


3.046 


-144.34 


5779.728 


5780.410 


DIB 


0.3405 






0.682 


-20.52 


5796.444 


5797.030 


DIB 


0.0440 






0.586 


-15.46 


6139.673 


6143.063 


Nel(l) 


0.0282 


-0.350 


16.61-18.62 


3.39 


-150.68 


6195.200 


6195.990 


DIB 


0.0756 






0.790 


-23.38 


6202.191 


6203.060 


DIB 


0.0966 






0.869 


-27.16 


6269.073 




UN 


0.0691 










6282.734 


6283.86"* 


DIB 


0.7319 






1.126 


-38.89 


6398.61 


6402.246 


Nel (1) 


0.0543 


0.360 


16.61-18.54 


3.636 


-155.48 


6612.857 


6613.630 


DIB 


0.1268 






0.773 


-20.19 


6637.522 


6641.031 


OII(4) 


0.0238 


-0.884 


23.40-25.27 


3.509 


-143.64 


6717.943 


6721.388 


011(4) 


0.0963 


-0.610 


23.43-25.27 


3.445 


-138.89 


7766.623 


7771.944 


OI(l) 


0.576 


0.320 


9.14-10.74 




blend 


7769.416 


7774.166 


OI(l) 


0.434 


0.170 


9.14-10.74 




blend 




+ 7775.388 


OI(l) 




-0.050 


9.14-10.74 






8592.524 




UN 


0.265 










8646.868 


8648.28 


DIB 


0.466 






1.412 


-34.11 



''DIB 6283. 86 A is blended with telluric absorption Unes in this region. 
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Table 3. Emission lines in IRAS 22023+5249 



/lobs. 


/llab. 


Went. 


w, 


log (gf) 


X 


lA^I 




(A ) 


(Aj 




(A ) 




(cV) 


(A) 


km 


4411.134 


4413.601" 


Fen(32) 


0.0600 


-3.870 


2.67-5.48 




blend 


4811.856 


4814.55 


[FeII](20F) 


0.0409 






2.694 


-153.00 


4902.59 


4905.35* 


[Fell](20F) 


0.0287 








blend 


5038.313 


5041.024 


Sill(5) 


0.1095 


0.291 


10.06-12.52 


2.711 


-146.47 


5053.306 


5055.984 


Sill(5) 


0.2758 


0.593 


10.07-12.52 




blend 




+ 5056.317 


SiII(5) 




-0.359 


10.07-12.52 






5088.76 




UN 


0.0244 










5155.834 


5158.81 


[FeIl](19F) 


0.0985 








blend 


5191.418 


5193.909 


Felll(5) 


0.0594 


-2.852 


8.65-11.04 




blend 




+ 5194.384 


Felll(5) 






8.65-11.04 






5195.026 


5197.929 


Fel(1091) 


0.1135 


-1.640 


4.30-6.68 


2.903 


-152.68 


5197.427 




UN 


0.084 










5232.982 




UN 


0.0412 










5240.659 


5243.306 


Felll(113) 


0.0997 


0.405 


18.26-20.62 




blend 




+ 5243.773 


Fel(1089) 




-1.150 


4.25-6.62 






5258.712 


5261.61 


[Fell](19F) 


0.0520 






2.898 


-150.36 


5270.327 


5273.38 


[Fell](18F) 


0.0481 






3.053 


-158.81 


5273.648 




UN 


0.0487 










5279.665 




UN 


0.0756 










5282.071 




UN 


0.0277 










5286.85 




UN 


0.0632 










5296.105 


5299.044 


01(26) 


0.1163 


-2.140 


10.98-13.32 




blend 


5297.325 




UN 


0.0366 










5299.853 




UN 


0.0368 










5751.43 


5754.8 


[NII](3F) 


0.0436 








weak 


5830.998 


5834.06 


Fell(165) 


0.0705 


-3.738 


5.57-7.69 


3.062 


-142.58 


5917.114 


5920.124 


Felll(115) 


0.1332 


-0.034 


18.78-20.87 


3.01 


-137.66 


5926.609 


5929.685 


Felll(114) 


0.0964 


0.351 


18.50-20.59 


3.076 


-140.75 


5950.424 


5953.613 


Felll(115) 


0.0989 


0.186 


18.78-20.86 


3.226 


-147.69 


5954.336 


5957.559 


SiII(4) 


0.1449 


-0.301 


10.06-12.14 


3.223 


-147.43 


5975.851 


5978.930 


SiII(4) 


0.4134 


0.004 


10.07-12.14 


3.079 


-139.62 


5996.369 


5999.70 


Alll(93) 


0.1843 




15.52-17.57 




blend 




+ 5999.83 


Alll(93) 






15.52-17.57 






6029.382 


6032.67 


Fel(1082) 


0.2320 




4.20-6.25 


3.288 


-148.64 



«FeII(32) 4413.601A is blended with 011(5) 4414.899A absorption feature. ''[FeII](20F) 4905.35A is blended with 011(28) 
4906.88A absorption feature. 



satellite also indicates considerable circumstellar reddening (E(B-V)ciicumsteiiar = 0.44, Gauba & 
Parthasarathy 2003). 
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3.3. Distance estimates 

Assuming 122023 to be a normal population I star, we may adopt the simple approximation given 
in Lang (1980) for the velocity structure of the Galactic disk which expresses the heliocentric radial 
velocity as a function of distance and Galactic coordinates of the star : 
Vo = -5vo + A D Sin 21 Cos^b [kms"i], 

where, v© is the heliocentric radial velocity of 122023; Svq is the correction term to convert from 
heUocentric radial velocity to velocity in the Local Standard of Rest (visr) and is 13.19 kms~^ for 
122023; Oort's constant, A ^ 15 km s~^ kpc"^; D is the distance to the star in kpc and the Galactic 
coordinates for 122023, (l;b) = (99.30; -1.96). Using Vq = - 148.31 km"' (Sec. 3.4), we derive a 
distance of 28.2 kpc to 122023. Such a large distance is unphysical and suggests that 122023 lies 
beyond the edge of our Galaxy. It is therefore, unlikely that 122023 is a normal population I B 
giant/supergiant star. 

On the post-AGB evolutionary tracks of Schonberner (1987), a post-AGB star with Tefj - 
24000K and log g = 3.0 (Sec. 3.8) has a core mass, = 0.565 Mq. If 122023 is a post-AGB star, 
then using the relation between core-mass and quiescent luminosity maximum (Wood & Zarro 
1981) for AGE stars, we estimate Mboi = -4.30. Applying the bolometric correction (= -2.43) for 
a Bim star (Lang 1992) we derive My = -1.87. The observed V magnitude (mv= 12.52; Hog et al. 
2000) of the star was corrected for extinction using E(B-V) = 0.95 (Ay = 2.94). Distance modulus 
method then yields a distance of 1.95 kpc to the star. 

3.4. Radial velocity 

Heliocentric radial velocities have been derived from wavelength shifts of the well defined absorp- 
tion and emission Unes (Tables 2, 3, 4 and 5). The mean heUocentric radial velocities from the 
absorption and emission lines (Tables 2 and 3) are -148.31 ± 0.60 kms"^ and -144.13 ± 0.72 
kms"' respectively. Radial velocity measurements of the forbidden Unes have been excluded in 
estimating the mean. The quoted errors refer to the probable errors of estimation. Fig.2 shows the 
radial velocity trends with respect to the equivalent widths (W^) and lower excitation potentials 
(LEP) of the absorption and emission lines respectively. 

The mean heUocentric radial velocity of the [N II], [O I] and [Fe II] lines is -152.90 ± 0.96 
kms~'. [SII]6717.0A and6731.3A lines in 122023 have a markedly difl'erent heUocentric velocity 
corresponding to -171.93 ± 1.36 kms"^ The difl'erent radial velocities argue for a non-spherical 
nebula. 

3.5. Wind velocities and mass loss rate from the P-Cygni profiles 

The BaUner lines, Hq,, H^g, Hy, a few He I, CII and Fe III lines in 122023 show P-Cygni behaviour 
(Fig. 1, Table 4) indicating ongoing mass-loss. Wind velocities were estimated from the blue ab- 
sorption edges of the well defined and unblended P-Cygni profiles (Table 4). The wind velocities 
of the Fe 111(5) Unes are markedly difl'erent from those of the other species. However, the wind 
velocities do not show any obvious trend with the lower excitation potentials or ionization poten- 
tials of the species. The mean wind velocity from the Balmer, He I and C II lines is 187.48 kms"' 
and that from the Fe 111(5) lines is 149.26 kms This difference in velocities indicates a deviation 
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Fig. 2. Radial velocity trends of the absorption and emission lines in IRAS22023-I-5249. Radial 
velocity measurements of the forbidden lines have not been plotted. 



from spherical synnmetry - possibly a bipolar morphology and the presence of a dense equatorial 
torus (e.g. Welch et al. 1999; Sahai et al. 2005). HST images may reveal the true morphology of 
the object. 

The Balmer Unes, Hq,, Hyg and Hy in 122023 are shown in Fig. 3 The velocities are in the 
hehocentric frame. Detailed modelhng of these hnes to derive the mass loss rate would be the 
subject of a future paper. The equivalent widths of the Hq, emission components are related to 
the mass loss rates in OB stars (Leitherer 1988). The Hq, emission component in 122023 has an 
equivalent width (W^ = 37.8A ) comparable to that of the 081 star, HD152408 (W^ = 34.7A ). 
From this, we estimate a mass loss rate of 1.23x10"^ Moyr"' (Leitherer 1988). However, this value 
may be treated with caution. In 122023, the large equivalent width of the emission component 
may be due to a large amount of gas in its circumstellar envelope or the presence of a bipolar 
envelope (Volk et al. 2004) and may not be directly related to the mass loss rate. 



3.6. Diffuse interstellar bands (DIBs) 

DlBs 5487.43A , 5780.410A , 5797.030A , 6195.990A , 6203.060A , 6283. 86A , 6613.630A and 
8648. 28 A were identified in the spectrum of the star (Table 2). They have also been identified in 
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Fig. 3. Observed Balmer profiles (dotted) in 122023 show P Cygni behaviour. The velocities are in 
the heliocentric frame. 



other hot post-AGB stars e.g. 1RAS01005+7910 (Klochkova et al. 2002), IRAS 13266-5551 & 
IRAS 173 1 1 -4924 (Sarkar et al. 2005). 

DlBs may have their origin in both the interstellar and the circumstellar media. A lower and 
upper limit on the predicted velocities of interstellar lines may be obtained by substituting D = 
and D = 1.95 kpc (the assumed distance to 122023) in the formula stated in Sec. 3.3 above (see e.g. 
Bakker et al. 1996). This yields a range of -13.19 kms"^ to -22.53 kms"^ This may however be 
treated as only a rough estimate. 

A comparison of the heliocentric radial velocities of the identified DIBs with the estimate above 
suggests that all of these fines most likely originate in the intersteUar medium. From the strength of 
the band at 5780.410A (W^ = 340.5 mA ), we estimate interstellar E(B-V) ^ 0.63 (Herbig 1993). 



3.7. Na I D2 and Na I Di lines 

Six components were identified in the Na I D2 and Na I Di lines (Fig. 4, Table 5). The velocities 
of absorption component 1 and emission component 2 are comparable with the mean hefiocentric 
radial velocities of the absorption and emission lines in the star (Sec. 3.4) suggesting that compo- 
nent 1 is of photospheric origin and component 2 arises in an extended envelope around the central 
star Comparing the heliocentric velocities of absorption components 3, 4 and 5 with those of DIBs 
observed in the spectrum of the star we may infer that these components originate in the interstel- 
lar medium. Component 6 is observed in emission with a velocity very different from the envelope 
velocity. This component may arise in a disk or in outflows around the central star. The velocity of 
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this component is comparable with the expansion velocity estimated from the nebular hnes (Sec. 
3.9.2). 

3.8. Atmospheric parameters and abundances from absorption iine spectrum 

The detection of He I and Si III absorption lines in addition to the C 111 lines and the absence of 
He 11 absorption lines, indicates a BO - B 1 supergiant spectral type for the central star. Comparison 
of the UV(IUE) spectrum of 122023 with standard stars suggests that the star is a B2-supergiant 
(Gauba & Parthasarathy 2003). Similar to our analysis of the high resolution optical spectrum of the 
hot post-AGB star, IRAS13266-5551 (Sarkar et al. 2005), we used Kurucz's WIDTH 9 program 
and the spectrum synthesis code, SYNSPEC (Hubeny et al. 1985) along with solar metalhcity 
Kurucz (1994) model atmospheres to derive the atmospheric parameters and elemental abundances 
of 122023 under the LTE approximation. 

The largest number of absorption lines in 122023 are those of O II and N II. We derived the 
oxygen and nitrogen abundance with the WIDTH 9 program for various combinations of Teff, log 
g and ft- We covered 18000 K < Teff < 25000K. From the Kurucz (1994) model atmospheres, the 
log g value was limited to a minimum of 3.0. For each combination of these parameters, we then 
synthesised the spectrum using SYNSPEC. The best fit to the observed spectrum was obtained for 
Teff - 24000 K, log g = 3.0, kms"'. Since strong lines are usually affected by microturbu- 

lence, the use of these lines in determining the atmospheric parameters of the star may contribute 
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to systematic errors. In estimating the atmospheric parameters and abundances, we excluded lines 
with > 200 mA . Line blends were also excluded from this analysis. 

3.8.1. He I lines 

Several absorption and P-Cygni type He I lines were identified in the spectrum of 122023. The 
absorption fines are blends with the exception of the strong 4387.929A He 1(51) (W^ = 372.2mA ) 
line and the He 1(50) 4437.551 A line (W^ = 78.6mA ). For the derived atmospheric parameters of 
the star, using spectrum synthesis, we estimated the helium abundance from the He 1(50) line. 

3.8.2. C II and C III lines 

Both C II and C III absorption lines were detected in the spectrum of the star. However, the number 
of these lines in 122023 is few. Further, all the identified carbon fines are blends with the exception 
of the CII(16) fine at 5151.085A . The carbon abundance was therefore estimated using spectrum 
synthesis alone. The region of the CII(16) lines (~ 5132A - 5151A ) was used for this purpose. 

3.8.3. on lines and O I triplet 

The largest number of absorption lines in the spectrum are those of O II. The O II lines with > 
200mA may be sensitive to NLTE effects. Such strong lines were neglected while estimating the 
oxygen abundance using WIDTH9. For the derived atmospheric parameters and oxygen abundance, 
using the SYNSPEC code, we could not obtain a good fit to such fines. 

The (total) equivalent width of the O I triplet in the spectrum of 122023 is 1.01 A . This is 
comparable to the 0.95 A equivalent width of the O I triplet in the spectrum of the B1.5Ia hot 
post-AGB star, LSn+34°26 (Garci'a-Lario et al. 1997b; Arkhipova et al. 2001). The O I triplet at 
/17773A is known to be sensitive to NLTE effects. Using log e(0) = 8.90 derived from the O II 
lines in the spectrum of 122023 (Table 6a), we could not obtain a fit to the O I triplet. 

3.8.4. N II and Ne I lines 

Several N II and two Ne I fines were identified in 122023. The abundances of these lines were 
estimated using WIDTH9. Once again, the mismatch between the strong (W^ > 200mA ) observed 
and synthetic N II lines indicates the need for a NLTE analysis of the spectrum of this star. 

3.8.5. Mg II lines 

Only one Mg II line could be identified in the spectrum of the star. This line is blended with 
Al 111(8). Since the Al IH abundance in 122023 is uncertain (see below), we did not attempt to 
estimate the magnesium abundance from the blended 448I.I26A Mg 11(4) fine. 

3.8.6. Al III lines 

Only four Al III lines could be identified. Three of these lines are blends with other species. We 
therefore, estimated the aluminium abundance from the single 5722.730 A Al 111(2) line using 
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spectram synthesis analysis and derived log e(Al)=6.79. This abundance from a single Une with 
W^=78.6mA may be treated as an upper limit. 

3.8.7. Si III lines 

The derived silicon abundance, [Si/HJ=-0.12 suggests that 122023 may be sUghtly metal deficient. 

3.8.8. Iron lines 

Iron abundance could not be estimated since the iron fines in 122023 appear only in emission or 
sfiow P-Cygni profiles. 

3.8.9. Error estimates 

The standard deviation (cr) which measures the scatter in the abundances due to individual fines 
of a species was estimated using WIDTH9 (Table 6a). The true error, a I y/n, would be smaller for 
species with a greater number of fines (n). 

Table 6b gives the uncertainties in the abundances due to uncertainties in the model atmospheric 
parameters, taking ATejf=-1000K, Alog g=+0.5 and A^t=+1 kms"^ The quadratic sum of the 
uncertainties due to the model parameters is given by cr^. 

3.9. Analysis of the emission line spectrum 

Several permitted and forbidden enussion fines were identified in the spectrum of the star (Table 
3). Nebular parameters and expansion velocities were determined using the forbidden lines. 

3.9.1 . Nebular parameters 

In the absence of a flux calibrated spectrum for 122023, it is not possible to obtain the absolute 
fluxes in the observed emission lines. However, reliable emission line flux ratios may be deduced 
by combining the observed equivalent widths (Table 3) with estimates of the stellar continuum 
flux distribution in the regions of the emission lines. The latter were obtained for the derived atmo- 
spheric parameters of the star (Sec. 3.8 ) using the SYNSEPC code and Kurucz model atmospheres. 
The emission line fluxes thus estimated are free from the effects of both intersteUar and circumstel- 
lar redderung. 

The [S II] /16717//16731 fine ratio is an electron density diagnostic and the [N II] 
(/16548+/16583)//15755 fine ratio is sensitive to electron temperature. In 122023, the [N II] 
6583.6A emission is blended with the emission component of the C 11(2) 6582.88A P-Cygru 
profile. However, comparing the 6582.88A C 11(2) profile with the 6578.05A C 11(2) P-Cygni 
profile (Fig. 1), we may conclude that the contribution of C 11(2) to the [N 11] emission profile is 
negligible. Using the NEBULAR analysis package under IRAF, we obtained Te Vs. Ng contours 
for the observed |S 11] and [N II] diagnostic ratios of 0.5 and 166.7 respectively. From the inter- 
section of the contours we obtained Te=7000 K and Ne=1.2xl0'* cm"-'. The high electron density 
is comparable to that measured in the young, compact PN, Hen3-1357 (Parthasarathy et al. 1993, 
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Bobrowsky et al. 1998) which evolved from the hot post-AGB stage into a PN in 20-30 yrs. 
(Parthasarathy et al. 1995). 

We could not derive the nebular C, N, O abundances which could then have been compared 
with the photospheric abundances to estimate the amount of material lost by the star during nebular 
formation and the chemical composition of the nebula. Such a calculation requires an estimate of 
the Hp emission line flux. However, Hp in 122023 shows a P-Cygni profile and it is not possible to 
estimate the nebular emission from this profile. 

3.9.2. Expansion velocities 

Expansion velocities were estimated from the FWHM of the unblended [Ol], [Nil] and [SII] lines 
using Vexp=0.50 FWHM (Table 7). This approximation is valid when emission is confined to a 
thin spherically symmetric shell. However, 122023 appears to have a bipolar morphology in recent 
ground based high resolution images (Volk et al. 2004). Further, the observed [OIJ 6300. 23A and 
6363. 88A Une profiles appear to be asymmetric. Though no obvious Une split is observed in the 
weak [OI] lines, their asymmetric nature may indicate the presence of a red and a blue component. 
This may explain the discrepancy between the expansion velocities estimated from [OI] and [Nil] 
and [SII] Unes in 122023, the former being nearly twice that of the latter two species. The mean 
nebular velocity based on the [Nil] and [SII] lines is 17.5 kms"' . 

4. Discussion and conclusions 

In this paper, we present the first high resolution optical spectrum and detailed line identification of 
the B-type star, LS III +52°24, identified as the optical counterpart of the hot post-AGB candidate 
122023. That 122023 is not a normal population I B star is also suggested by the large heUocentric 
radial velocity of the star (-148.31 ± 0.60 kms"^) measured from the absorption Unes in its high 
resolution optical spectrum. It is more Ukely that 122023 is a post-AGB star belonging to the old 
disk population. The presence of absorption lines of He I, C III, Si III alongwith [N II], [O I] and 
[S II] emission lines indicate a low excitation nebula surrounding an early B-type star. P-Cygni 
lines of hydrogen. He I, C II and Fe III provide evidence for ongoing post-AGB mass loss. 

As a first approximation,using LTE analysis, we estimated Teff - 24000 K ± 1000 K, log g - 
3.0 ± 0.5 and = 7 ± 1 kms"'. The O 1 triplet and strong He 1, O 11 and some N 11 absorption lines 
indicate that NLTE effects are significant in the star. The derived CNO abundances were compared 
with average CNO abundances for main- sequence B- stars from the Ori OBI association (Table 
6a; KiUan 1992). The CNO abundances suggest that 122023 is an evolved star. We estimated C/O 
~ 0.48. A spectrum with better S/N ratio and NLTE analysis of the absorption lines may provide 
improved estimates of the atmospheric parameters and abundances. Besides high resolution UV 
spectra of the object may reveal a greater number of absorption fines in the UV which would help 
in a more accurate determination of these parameters. 

Even though the IRAS cool and hot post-AGB stars show supergiant like spectra indicating an 
evolutionary sequence in the transition region from the AGB to the PN stage, they seem to show 
significant differences in chemical composition. In the high galactic latitude hot post-AGB stars 
severe carbon deficiency is detected indicating that they left the AGB before the third dredge-up 
(Conlon et al. 1991; McCausland et al. 1992; Moehler & Heber 1998). Similar carbon deficiency 
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is also detected in the hot post-AGB stars in globular clusters (Moehler et al. 1998; Mooney et al. 
2004; Jasniewicz et al. 2004; Thompson et al. 2006). The only exception to this pattern is the field 
hot post-AGB star, IRAS 01005+7910 which showed an overabundance of carbon (Klochkova et 
al. 2002). In contrast, amongst IRAS selected cool post-AGB stars, not a single post-AGB star 
is detected with severe carbon deficiency. The F-type post-AGB stars with 21 -micron emission 
show an overabundance of carbon and s-process elements indicating that they have experienced 
third dredge-up and evolved from the carbon star stage on the AGE. 

Evidence for bipolar nebulosity was seen in ground based H2 images of 122023 (Volk et al. 
2004). Bipolar proto-planetary nebulae (PPNe) are usually attributed to binary systems. In our 
high resolution optical spectrum, we did not find any evidence for binarity of the hot central star. 
HST WFPC2 images may provide better insight into the morphology and evolutionary nature of 
this object. 

Far-IR (IRAS) colours similar to PPNe, H2 emission arising close to the central star with 
possible bipolar morphology (Volk et al. 2004), large heliocentric radial velocity of - 148.3 1 kms"^ 
and absence of large amplitude irregular light variations similar to that observed in luminous blue 
variables (LBVs) indicate that 122023 cannot be classified as a LBV candidate. 122023 may show 
low amplitude irregular light variations similar to that observed in SAO 85766, LSII +34°26 and 
other hot post-AGB stars (Arkhipova et al. 2006). 
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Table 3. Emission lines in IRAS 22023+5249 
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6300.23 


[0I](1F) 


0.2832 






3.461 


-149.93 


6337.276 


6340.58 


NII(46) 


0.0479 


-0.192 


23.23-25.18 


3.304 


-141.46 


6343.686 


6346.86 


Nll(46) 


0.5132 


-0.901 


23.22-25.18 




blend 




+ 6347.109 


Sill(2) 




0.297 


8.12-10.07 






6353.584 


6357.0 


NII(46) 


0.0511 




23.23-25.18 


3.416 


-146.34 


6360.279 


6363.88 


[0I](1F) 


0.0677 






3.601 


-154.89 


6367.944 


6371.371 


Sin(2) 


0.1989 


-0.003 


8.12-10.06 


3.427 


-146.49 


6458.428 




UN 


0.1059 










6544.492 


6548.1 


[N1I](1F) 


2.817 






3.608 


-150.43 


6579.838 


6583.6^^ 


[N1I](1F) 


8.534 








blend 


6712.772 


6717.0 


[S1I](2F) 


0.5017 






4.228 


-173.96 


6727.154 


6731.3 


[SII](2F) 


1.007 






4.146 


-169.91 


6848.05 


6851.634 


Fel(34) 


0.1321 


-5.320 


1.61-3.41 


3.584 


-142.06 


6998.28 


7001.93 


01(21) 


0.0786 




10.94-12.70 




blend 




+ 7002.22 


01(21) 






10.94-12.70 








7231.330'' 


Cll(3) 




0.043 


16.32-18.03 




blend 




7236.420'' 


CII(3) 




0.299 


16.32-18.03 




blend 




7316.282'' 


atmos. 












7373.821 




UN 


0.1893 










7458.676 




UN 


0.1119 










7462.426 




UN 


0.1229 










7464.227 




UN 


0.122 










7892.292 




UN 


0.134 










7993.283 


7993.332 


atmos. 


0.0534 










8218.700 


8223.128'' 


NI(2) 




-0.390 


10.32-11.83 




blend 


8237.897 


8242.389'' 


NI(2) 




-0.380 


10.33-11.83 




blend 


8283.265 




UN 










blend 


8441.864 


8446.359 


01(4) 


3.097 


0.170 


9.51-10.98 




blend 




+ 8446.758 


01(4) 




-0.050 


9.51-10.98 






8745.796'' 




UN 












8857.744'' 




UN 













^[NII](1F) 6583.6A is blended with the emission component of CII(2) 6582.88A P-Cygni profile. ''These 
emission lines are weak and are blended with the atmospheric absorption lines in this region. 
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Table 4. P-Cygni lines in IRAS22023+5249. Equivalent widths of the absorption and emission 
components of the P-Cygni profiles are given. Wind velocities are estimated from the blue absorp- 
tion edges of the P-Cygni profiles. 





Ident. 


(absorption) 


(emission) 


log (gf) 


y 


Wind Velocity 


(A) 




(A) 


(A) 




(eV) 


km s"' 




XJ 




101^ 




1 n 1 1 oj. 


IcS/.oi 








n nQQ7 

u.uyo / 


7 7 1 S 
— Z.Zl o 




blend 


AA^ 1 n 1 Q 


V(^^'\^(A\ 




u.ujud 


7 ^^77 
— Z.J / z 


8 9zl 11 OA 


weak 


zlA71 All 


\^p^^(\A\ 

riCi.yi'-r ) 


U.J J)70 


U.UoZo 




70 Q'^ 9^ 79 
ZU.7J— Zj. /Z 


blend 


A.AA1\ fStl 


\Af^^(^A\ 






n 808 


90 0^ 9'^ 79 








U. 


f% IIQ 


—0 090 


10 10—19 Id 
L\j, ly — iz. /*+ 


— 109 fi"? 


JUi J.U / o 


nciv^H- ^ 




U. JUH- / 


—0 890 


90 ^^0— 9^ 07 


— 18'? ^0 

— lOJ. J7 






n 1 ^Ai 


n nA7^^ 


— Z.J J / 


a 1 1 00 

O.OJ— 1 1 .yJy 


weak 


'^nSA 701 




U.U J 1 y 




— Z. J7U 


a AS 1 1 00 

o.Oj— i 1 .\Jy 


weak 


'^1 77 ^87 




n 1 OAS 


n 1 f\^A 


7 718 


a 1 1 07 

o.Dj— i l.U / 


1 AQ 


CI 111 

J 1 JD. Ill 


FcliH^J ) 


U.i lUJ 




9 01 8 


8 AJ. 1 1 04 


ld.0 Id. 


Jo / j.Dio 


r\S-\\\- 1 ) 


U.JJUU 


1 8700 
l.o /yu 




90 87 99 07 
ZU.o /— ZZ.7 / 


blend 




Hel(ll) 








90 87-99 07 




+ 5875.989 


Hel(ll) 








20.87-22.97 




6562.797^ 


H„ 


0.0443 


37.86 


0.710 


10.19-12.08 


-185.36 


6578.050 


CII(2) 


0.1977 


0.2779 


-0.026 


14.44-16.32 


-191.55 


6582.880^ 


cn(2) 






-0.327 


14.43-16.32 


blend 


6678.154 


Hel(46) 


0.6693 


0.8538 


0.329 


21.20-23.06 


-182.16 


7065.176 


Hel(lO) 


0.3406 


1.3600 


-0.460 


20.95-22.70 


blend 


+ 7065.707 


Hel(lO) 






-1.160 


20.95-22.70 





"The absorption component of Felll(l 14) 4419.596A is blended with 011(5) 4416.975 absorption feature. ''The 
emission components of the Hp and H^, profiles have broad wings. Gaussian fits to the absorption and emission 
components of these profiles could not be obtained. Using IRAF, the equivalent widths of these components 
were estimated by subtracting the linear continuum between the points of interest and summing the pixels with 
partial pixels at the ends. ^ CII(2) 6582.88A P-Cygni profile is blended with [Nn](lF) 6583.6A . 



Table 5. Absorption (a) and emission (e) components of Na I D2 (5889.953 A ) and Na I Di 
(5895.923 A ) fines in the spectrum of IRAS22023+5249 (LSIII +5224). are the equivalent 
widths of the components and Vr are the respective heliocentric radial velocities. 







IRAS22023+5249 






NalDz 


NalDi 


Component 




/lobs. 


w. 


V, 


^obs. 




V,- 






(A) 


(A) 


(km s-i) 


(A) 


(A) 


(km s-') 


1. 


a 


5886.341 


0.0559 


-169.10 


5892.323 


0.0529 


-168.31 


2. 


e 


5886.811 


0.1314 


-145.17 


5892.763 


0.0277 


-145.92 


3. 


a 


5888.530 


0.1395 


-57.61 


5894.500 


0.1355 


-57.54 


4. 


a 


5888.903 


0.5054 


-38.61 


5894.881 


0.3684 


-38.15 


5. 


a 


5889.500 


0.6384 


-8.20 


5895.455 


0.6241 


-8.94 


6. 


e 


5889.918 


0.1221 


+ 13.09 


5895.902 


0.1471 


+13.8 
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Table 6a. Derived abundances for 122023. The values are for log e(H) = 12.0. n refers to the number 
of Unes of each species used for the abundance determination. For comparison we have Usted the 
solar abundances (log e(X)o)and average abundances for main sequence B-stars from the Ori OB 1 
association (KiUan, 1992). 







IRAS22023+5249 




J-VJ.Cll.ll lSC'UUVIIV'C' 




(Teff 


=24000 K, log g=3.0, ^t=7 km s"') 




B-stars, Ori OBI 


X 


n 


log e{X) 


cr 


[X/H] 


loge(X)o 


log e(X) 


He I 


1 


11.04+ 




+0.11 


10.93 


11.04 


CII 


4 


8.58+ 




+0.19 


8.39 


8.23 


Nn 


7 


8.36 


0.33 


+0.44 


7.92 


7.72 


on 


18 


8.90 


0.42 


+0.21 


8.69 


8.60 


Nel 


2 


9.02 


0.14 


+0.94 


8.08 




Aim 


1 


< 6.79"'" 






6.47 




Sim 


3 


7.43 


0.59 


-0.12 


7.55 





The abundances were derived using Kurucz's 'WIDTH9 program. 

+These values were derived from spectrum synthesis analysis using the SYNSPEC code. 



Table 6b. Uncertainties in the derived abundances, Alog e(X), due to uncertainties in the model 
atmospheric parameters 



Element 


ATeff 


Alogg 


A^t 






-lOOOK 


+0.5 


+1 km s"' 




N 


-0.14 


-0.18 


-0.05 


0.23 





+0.05 


+0.08 


-0.08 


0.12 


Ne 


-0.15 


-0.26 


-0.01 


0.30 


Si 


-0.04 


-0.06 


-0.07 


0.10 



Table 7. Expansion velocities 



Went. 


-llab 


FWHM 


V 

^ exp 




A 


A 


km s"' 


6300.23 


[01](1F) 


1.334 


31.76 


6363.88 


[01](1F) 


1.470 


34.65 


6548.1 


[N11](1F) 


0.844 


19.34 


6717.0 


[SII](2F) 


0.747 


16.68 


6731.3 


|S1IJ(21'") 


0.739 


16.47 



[NII](1F) 6583. 6A line is blended with the emission component of CII(2) 6582.88A P-Cygni profile and has 
not been used to estimate the expansion velocity. 




Fig. .1. Optical spectram of IRAS22023+5249 (LS III +52°24) 
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Fig. 1. Optical spectrum of IRAS22023+5249 (LS III +52°24) contd.. 
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Fig. 1. Optical spectrum of IRAS22023+5249 (LS III +52°24) contd.. 
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